Na 4 Ir 3 08 as a 3D spin liquid with fermionic spinons 
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Spin liquid states for spin-1/2 antiferromagnetic Heisenberg model on a hyperkagome lattice 
are studied. We classify and study flux states according to symmetries. Applying this model to 
Na4lr3C>8, we propose that the high temperature state may be described by a spinon Fermi surface, 
which forms a paired state with line nodes below 20 K. The possible mixed spin singlet and spin 
triplet pairing states are discussed according to the lattice symmetry which breaks inversion. 
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A spin liquid is a spin system where quantum fluctu- 
ations dominate its low energy behavior, thereby a long 
range spin order can not be established even at zero tem- 
perature. A spin liquid is expected to have exotic proper- 
ties such as spinons carrying S =1/2 excitations. After a 
long search, promising examples in two dimensions (2D) 
have been identified. [lj Recently a spinel related oxide, 
Na4lr 3 0g, was proposed as the first candidate for a 3D 
spin liquid Temperature dependent spin susceptibil- 
ity around room temperature yields an antiferromagnetic 
(AFM) Curie-Weiss constant &w ~ 650 K, and there is 
no anomaly indicative of long range spin ordering down 
to 2K. Some kind of phase transition or cross-over seems 
to occur at a temperature T c ~ 20K in that the specific 
heat Cy divided by T shows a rather sharp peak. On 
the other hand, the spin susceptibility x {T) is almost 
temperature independent. Using the experimental value 
of spin susceptibility \ and specific heat ratio 7 at the 
specific heat peak at ~ 20K, we find that the Wilson 



ratio Rw 



g k BX 



of the material is 0.88, which is very 



close to that of a Fermi gas where Rw is unity. There- 
fore, for a wide range of temperature T c < T < 6\y, the 
system seems to be a Fermi liquid of spinons. Below T c 
the specific heat decreases to zero as Cv ~ T 2 , indicating 
a line nodal gap in the low lying quasiparticle spectrum. 
However, this picture needs to be reconciled with the ob- 
servation that the spin susceptibility \ remains almost 
constant. In this letter, we will explain such phenomena 
based on a fermionic spin liquid picture. 

The spins come from low spin 5d 5 Ir 4+ ions and form 
a 3D network of cornered shared triangles. This gives 
rise to a S = 1/2 AFM coupled spin system on a hyper- 
kagome lattice (see FigJT}. It should be cautioned that 
because of the large atomic number, spin-orbit coupling 
in Ir is expected to be strong. This question was recently 
addressed by Chen and Balents.Q We shall first discuss 
the case when spin-orbit coupling is negligible and com- 
ment on the strong coupling case later. 



A hyperkagome lattice is a cubic lattice with 12 sites in 
a unit cell. There are two types of hyperkagome lattices 
corresponding to space group PAi?>2 and P4 3 32 respec- 
tively. They differ just by chirality. [2| In this letter, we 
will focus on P4i32. The chiral nature of the lattice is 
clearly seen as follows. Each lattice site is connected to 
two bonds which forms a straight line (shown in green 
in Fig. 1(a)). Around each triangle, these straight lines 
define a chiral pinwhccl. A set of arrows on each bond 
of the triangle can now be uniquely defined as shown in 
Fig. 1(a). The space group P4i32 contains 24 symme- 
try operations; its corresponding point group is octahe- 
dral group 0[4|]. It is readily seen that all the hyperk- 
agome sites and bonds between nearest neighbor sites are 
equivalent Since a unique oriented triangle is attached 
to each bond, all oriented triangles are also equivalent. 

We start with a nearest-neighbor (NN) Heisenberg 
model with AFM exchange J and consider spin liquid 
states on such a hyperkagome lattice. The spin liquid 
state is described by a fermionic trial wave function [j| 
which is the ground state of a trial Hamiltonian, 



H, 



trial 



(i) 



where Hq is a tight binding Hamiltonian of spinons and 
Hpair describes pairing between two spinons with 

H = - ^,^4^°^ + h.c, (2) 

(%H,jv),a 

Hpair = -J2 A a p{k)cl a cl kf3 + h.c., (3) 

ka/3 

where i denotes unit cells, fi labels sites in a unit 
cell, (in,ju) denotes a pair of NN sites, c] (c^) de- 
notes a fermionic spinon creation (annihilation) operator, 
A a( a(k) is a gap function where a =|, | are spin indices. 
The spin liquid state is formed by Gutzwiller projection 
of the trial wavefunction into a state with no double occu- 
pancy, tifiju and other parameters in the trial Hamilto- 
nian should be viewed as variational parameters obtained 
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FIG. 1: (Color online) Hyperkagome lattices with space group 
P4i32. (a) 8 triangles giving rise to 3-site directed loops. 
Straight green lines indicate a chiral pin-wheel, (b) A 10-site 
directed loop. Eleven other 10-site directed loops can be gen- 
erated from it by symmetry operations, (c) A 16-site directed 
loop shown by red arrows. The green arrows show two ori- 
ented 10-site loops that are related by symmetry operations. 
The figure shows how to decomposite the 16-site loop into 
sums of two 10-site loops and four 3-site loops, (d) An il- 
lustration of Eq.Q, where /i(i) denotes a site in unit cell 
i. We unfold the lattice onto 2D and keep the topology, i.e. 
the surface enclosed by the third 10-site loop which makes up 
the rest of a closed surface is mapped to the rest of the 2D 
plane outside the figure. Arrows indicate the orientation of 
three 10-site loops related by symmetry. The left 10-site loop 
is shown in 1(b). 



by minimizing the ground state energy through a micro- 
scopic spin model. In this letter we treat the problem at 
the mean field level and consider spin liquid states with 
full P4i32 point group and translational symmetries. 

First we consider Hq- The hopping integral t^jv can 
be written as t^ju = \t\ e ~ lA ^- iv with = -Aj Vtitl . 

We shall identify all possible flux states with translational 
and lattice symmetries. 

All the loops on a hyperkagome lattice can be decom- 
posed into two kinds of "elementary" loops, the 3-site 
and 10-site loops (see FigO}. As explained earlier, the 
3-site loops are equivalent under symmetry operations 
with the specified loop directions (see fig.(l)). Similarly 
with the 10-site loops, once the direction of one loop is 
fixed, symmetry operators will generate all 12 loops with 
specified directions. Therefore all the 3-site and 10-site 
loops have same flux $3 and $10, respectively, for states 
that respect translational and P4i32 symmetry. 



It is also helpful to define the 16-site loop shown in 
Fig. 1(c). This loop forms a surface and repetition of this 
surface cuts through the whole lattice in the xz plane. All 
the 16-site loops form surfaces covering different lattice 
planes and can be transformed to each other by sym- 
metry operation. The 16-site loops can be decomposed 
into sums of four 3-site loops and two 10-site loops as 
shown in Figfljc). Since the green arrows of one of the 
10-site loops are opposed to the red arrows, the flux on 
the two 10-site loops cancel. Similarly, the <&3's also can- 
cel, so that $16 = 0. We conclude that flux arrangements 
which satisfy the lattice symmetry do not permit a net 
flux through the lattice. 

Next we consider a volume enclosed by three 10-site 
loops and two 3-site loops as shown in Fig. 1(d) projected 
to the 2D plane. Since the total flux through the closed 
surface must be integer multiples of 27T, we find 



2$ 3 = 3$io(mod27r). 



(4) 



The states with zero and 7r flux through the triangles are 
equivalent after projection. Moreover, the states with $3 
and 7r — $3 are the same state after projection. This is 
because if we make a particle hole transformation /ij = 
/j, we have — > — t*j, but the projection of the particle 
and hole states are the same at half- filling. [|| Thus we 
can restrict to —ir/2 < $3 < ir/2. There is only one time 
reversal invariant state $3 = $10 = which does not 
break lattice symmetry. At half filling, Fermi surfaces of 
this state appear in only three bands as shown in Fig. [5] 
The other bands are either fully filled or empty. 



<I> 3 =0. 6th band 



*,=0, 7th band 




FIG. 2: Fermi surfaces of the zero flux state at half filling. 

Breaking time reversal symmetry will lead to flux 
states with arbitrary fluxes $3 and $10 satisfying Eq.(j4]), 
and still requires <&i6 = 0. All of these states, with or 
without time reversal symmetry, have finite Fermi sur- 
faces at half filling (see Figs. [2] and [3]). 

To determine what is the most stable state, we perform 
a mean field theory calculation on the AFM Heisenberg 



Sj V with trial Hamiltonian 



model, H = JJ2 

Htriai = Ho, first neglecting the pairing terms. By stan- 
dard U (1) slave boson mean field theory, we have the 



relation t 



J 



ergy is Eo/Nt, = — 8t 2 /3J, where Nb is the number of 
bonds. We find that Eq reaches its minimum at $3 = 0. 



and the ground state en- 
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=51/4, 6th band 




FIG. 3: Fermi surfaces of flux states with broken time reversal 
symmetry at half filling. 



Phcnomcnologically we expect that Hq determines the 
physics of the spin liquid state at T > T c and a spinon 
pairing gap characterized by H pa i r is opened up at T < 
T c . The power law behavior Cy oc T 2 observed at T < T c 
indicates that the gap has line nodes on Fermi surfaces. 
For triplet pairing, group theoretical analysis tells us that 
a spin triplet pairing state on a cubic lattice can create 
only full or point nodal gapsQ- This would seem to im- 
ply singlet pairing. However, because of broken inversion 
symmetry on a hyperkagome lattice^, spin singlet and 
triplet pairing states can now mix together in the pres- 
ence of spin orbit coupling^. In terms of d- vector the 
gap function A Q , j a(k) can be written in matrix formQ, 



A(k) = i (d (k) a + d (k) ■ a) a y , 
with the energy dispersion of quasiparticles 



parity (spin singlet) and even parity (spin triplet) pair- 
ing states within an irreducible representation T of oc- 
tahedral group O. For simplicity, wc shall only consider 
s-wave, (or extended s-wave), p-wave and d-wave states. 
In this case the available possibilities are Ti (s+p-wave), 
T3 (d + p-wave) and (d + p-wave) . We shall also limit 
ourselves to spin rotation invariant states with (S) = 
for simplicity. Therefore d (k) x d* (k) = and d (k) is 
real apart from a overall phase factor. The lower energy 
gap is then = y/dod^ + d • d* {d^d + dod*\ and re- 
sults in a nodal surface (Ag = 0) in fc-space because of 
mixing between different parity states as we shall see in 
the following. A nodal line of gap in quasi-particle ex- 
citation spectrum is produced by the cross line between 
nodal surface and Fermi surface (£g = 0). 

First we consider ID representation Ti. In this case the 
spin singlet part do may be s-wave or extended s-wave, 
with a p-wave state for spin triplet part d. For s-wave, 
we have do(k) = A s and d(k) = A t (smk x x + sin k y y + 
sm k z z). A g = implies sin 2 fc^+sin 2 fcj,+sin 2 k z = \^-\ 2 
which produces a nodal surface. Notice that the nodal 
surface shrinks to a point when A s — > 0. For extended 
s-wave, we have do (k) = A s (cos k x + cos k y + cos k z ) and 
d(k) = A^sinfc^a; + sin k y y + sin k z z). The node surface 

l I A = 1 2 sin 2 fc^+sin 2 /c^+sin 2 k z 

is given by ^ r = -, r^f j—^ — rnf-. 

° J 1 At 1 (cos /c^+cos fcy+cos k z ) z 

Similar analysis can be applied to 2D representa- 
tion F3 and 3D representation T§. For representa- 
tion T3 (d + p-wave), we have do(k) = A s [2uiCosfc z — 
(ui + V3u 2 )cosk x - (ui - \fiu 2 ) cos k y ] and d(k) — 
At [2vi sin k z x — (vi + V3v 2 ) sin k x y — (vi — V3v 2 ) sin k y z] , 



(5) where u 2 



1 with iti(2) and v i{2) being 



real. For representation T5 (d+p-wave), we have do(k) 
A s [u x sin k y sin k z + u y sin k z sin k x + u z sin k x sin k y ) and 
- d(k) = A t [v x (smk z y + sink y z) + v y (sink y z + sinfc z x) + 

d d* + d • d* ± \/\dxd*\ 2 + (d*d + d d*f, v * ( sin M + sin k *v)\ > where u l + u l + u l = v l + v l + v l = 

1, with u x ( y z ) and v x t VtX \ being real. For all these states 

it is possible to have nodal surfaces which will shrink to 

points as A s — > 0. Thus we conclude that the pairing 

must be singlet or singlet with triplet admixture due to 

spin-orbit coupling in order to have node lines. 

Next we consider the spin susceptibility of such mixed 
states. We neglect spin-orbit interaction and consider S 
as 'real' spins at first. Singlet coupling predicts the spin 
susceptibility will go towards zero below T c , which is not 
observed experimentally. If both singlet and triplet are 
present and spin-orbit scattering is weak (in the sense 
to be defined later), we find the ^-dependent electronic 
contribution to spin susceptibility is given by *^ ^ = 
d Q d 'o+ d * d ' , d o d o+ d :^ Y (k- T) where i = x v z and 

d d*+d-d* T d d*+d-d* 1 1 /' wiieie t s £/! z > duu 

Y(k;T) is the fc-dependent Yosida function[9(. Assum- 
ing that the d vector is pinned to the lattice, for a 
polycrystalline sample we should average over all spa- 



. (6) 

where £k = £k — £k is the energy dispersion of spinons 
calculated from Hq and fi is the chemical potential. We 
have assumed £_k = £k in deriving E^, which is valid at 
half filling in presence of time-reversal symmetry. The 
unitary condition[9( is given by id x d* +d^d + d a d* = 0. 
It is obvious that the mixing of singlet and triplet states 
leads to nonunitary states even if the d-vector is real. 
For a nonunitary state, there are two branches of E^ 
with the energy gap determined by the lower branch. We 
will see how such mixing may lead to nodal line states 
later. In such pairing states, the U(l) gauge symmetry is 



t J 

icr^jcr' I 



However, a local 



broken by the spinon pairing (c, 

Z 2 symmetry, where Ci a ±Ci a together with a sign 
change in the gauge variable on the link remains intact, 
so that they can be classified as Z 2 spin liquids, fiol] 

Spin pairing states on a cubic lattice with inversion 
symmetry are classified in Ref.Q- What we need to do 
is to mix the pairing states as a linear combination of odd 



1 



tial directions, resulting in 



\do\ d 



3 |do| 2 + |d|- 



4 



| ffl s )y(T), where Y (T) is the Yosida function, Xs 
is the spin susceptibility below T c and xjv is the Pauli 
spin susceptibility at the normal state. Therefore ^ re- 



duces to | 



1?- 



+ |d| 



? at zero temperature. If the spin 



|; if the spin singlet 
Neither of these cases are 



triplet pairing dominates, ^ 
pairing dominates, — * 0. 
observed in experiment, where x changes little below T c . 

In order to explain the absence of change in x below 
20 K, we appeal to the effect of spin-orbit coupling. It 
is well known that in conventional BCS singlet super- 
conductors, the Knight shift hardly changes below T c for 
heavy elements such as Sn and Hg. ll| It was quickly 
realized that this is due to the destruction of spin con- 
servation due to the spin-orbit coupling. A clear expla- 
nation was given by Anderson [13], who introduced the 
notion of time reversed pairing states. Let us consider 
the imaginary part of the spin response function x"(<7> w )- 
If total spin is conserved, the dynamics is diffusive and 
x"{q,to) will have a central peak in oj space with width 
Dq 2 which goes to zero as q — > 0. Superconductivity 
gaps out all low frequency excitations, thus removing this 
central peak. By Kramers-Krong relation the real part 
x'(q = 0,ui — 0) vanishes in the superconducting ground 
state. In the presence of spin-orbit coupling, the total 
spin is not conserved, but decays with a lifetime t s . In 
this case x"{0. = 0) u ) nas a central peak with a width i. 
The superconducting gap A cuts a hole in the x"(w) for 
u> < A, but leaves the region lu 3> A intact, in agreement 
with physical expectation that the high frequency region 
should be unaffected by pairing. By Kramer-Kronig re- 
lation, x' wm be reduced, but if the spin-orbit coupling 
is sufficiently strong such that 



1 



» A 



(7) 



the reduction will be small, i.e. = 1 — 0(At s ). Equa- 
tion ([7]) is the criterion needed to have very little charge 
in the spin susceptibility below T c . 

Let us now discuss to what extent our model may be 
applicable to Na^^Os- The t^ g levels are split by crystal 
fields due to the presence of Na by an amount of order 
E%. Chen and Balents [3| made the important observa- 
tion that the tig levels may be considered as a L = 1 
multiplet and distinguish between the strong and weak 
coupling limits, depending on the spin-orbit energy A. In 
the strong coupling case 



A>£ 3 



(8) 



the spin and orbital degrees of freedom lock to form 
J = L + S and it turns out that J = In this case, 
the direct exchange between the Ir ions gives an isotropic 
Heisenberg model J Si ■ Sj, if S is interpreted as J. The 
g factor is 2 while superexchange via oxygen gives rise to 
anisotropy and Dzyaloshinski-Moriya (DM) terms. Since 



the experimentally measured g coupling is close to 2, 
Chen and Balents concluded that the Na^^Os system 
must be either close to weak coupling or strong coupling 
in the sense of Eq.®, and strong coupling seems more 
likely, given the size of A. Provided that the exchange 
path is dominated by direct exchange, our starting point 
with the isotropic Heisenberg model and the fermion rep- 
resentation model remains valid, provided that the spin 
labels in Eq.(2) are to be understood as J z — ±4. The 
correction terms are anisotropy and DM terms which are 
expected to be of order (E3/X) 2 w (g — 2) 2 in this case. 
These correction terms spoil the conservation of total J, 
giving rise to a decay rate analogous to — . In the limit 
At, <C 1 we again expect that the spin susceptibility 
is not much affected by the onset of pairing. This may 
be the case closest to experiment. We emphasize that 
the criterion for strong or weak spin-orbit coupling given 
by Eq.® is different from Eq.© which determines the 
magnitude of the drop in x- 

In summary, we studied fermionic spin liquid states 
described by trial Hamiltonian ((T|) on a 3D hyperkagome 
lattice, and classified all possible flux states, which all 
have Fermi surfaces. The broken inversion symmetry 
may lead to the mixing of spin singlet and triplet states 
in the spinon paired states, resulting in formation of 
line nodal gaps, which provides a natural explanation to 
the observed low temperature specific heat CV ~ T 2 in 
Na^^Os- The strong spin-orbit coupling explains to the 
observed spin susceptibility, which is more or less tem- 
perature independent at temperatures below ~ T c . The 
transition at 20 K is then interpreted as a transition be- 
tween a U(l) spin liquid to a Z2 spin liquid where spinons 
are paired at low temperature. This theory also predicts 
the appearance of superconductivity if the Na^^Os sys- 
tem can be doped. 
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